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PREFACE 

. Emphasis continues to bo placed on thf, uso of digital computers in 
solving nonlinear hydrodynamic and thermodynamic equations of fluid 
flow. This publication of the thermal pollution group at the University 
of Miami presents the solution of one such problem. This problem deals 
with the use of a numerical one-dimensional model in predicting the tem- 
perature profiles of a deep body of water. Although this model can be 
applied to most lakes, a s^acfflc site (Lake Keowee, S. C.) application 
has been chosen and described in detail. The programs are written in 
fortran V and could be modified by the user. Some of these modification?; 
are suggested either in the text or in the specific programs. 

A detailed derivation of the equations integrated has been left out; 
however, to improve readability of the final equations, the meaning of 
the terms and variables occurring in these equations are included. 

This research was performed at the thermal pollution laboratory at 
the University of Miami. Funding was provided by the National Aeronau- 
tics and Space Administration (NASA-KSC) and the Environmental Pro- 
tection Agency (EPA-RTP) . 
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ABSTRACT 



A user';? manual for a one- dimensional thermal modal Is described. 
The modal is essentially a set of partial' differential equations which are 
solved by finite difference methods using a high speed digital computer. 
The main equations integrated are discussed. The programs are written 
in fortran V and an example problem is discussed in detail. 
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SECTION 1 
INTRODUCTION 

It is important that the thermal behavior of heated discharges and 
their receiving basins be clearly understood. 

A numerical model that can be used for predicting the seasonal thermo- 
cline of a deep body of water Is very useful in studying the environmental 
impact of thermal discharges from power plants. This is not only required 
for existing power plants but also for planned units. Thus, a predictive 
capability is essential to the licensing procedure. Monitoring programs 
cannot sa'tlsfy these needs, but from time to time, play a vital role in the 
calibration and verification of mathematical models. 

* 

The one-dimensional, thermal numerical model, described in this 
manual, features the effects of area change with depth, nonlinear inter- 
action of wind-generated turbulence and buoyancy, absorption of radiative 
heat flux below the surface* thermal discharges and the effects of vertical 
convection caused by discharge. The main assumption in the formulation 
of this model is horizontal homogeneity. 

This model can be applied to most stratified deep bodies of water. 
This stratification has a seasonal cycle and is an important natural 
characteristic of a body of water. The body of water could bf divided 
into any number of slices. The temperature of each slice is predicted 
by the model. The surface slice exchanges heat with the environment 
of known climatic conditions wNle the bottom slice is assumed perfectly 
insulated. Condenser cooling water is extracted from any one of the 
slices and heated by the power plant. The discharge is injected into a 
slice of the same temperature as the discnage. 

The main function of the model is the prediction of the temperature 
profiles in a deep body of water for any number of annual cycles. How- 
ever, predictions connot be made on hourly basis - a feature usually 
handled by a more sensitive three-dimensional model. This is the main 
limitation of the model. 

The procedure used in writing this manual is as follows; 

Description and flow chart of the main program are given in Section 
3, where the subroutines are also described. In the next section, a 
list of the variables and dimensions are given. The next three sections 



•how how a typical run la prepared, executed and plotted. An axampU 
cm« is discuised in Appendix A, while Appendix B gives the fortran 
source program listings. 



SECTION 2 
RECOMMENDATIONS 



The main disadvantage of a one-dlmenslonal thermal modal Has In 
tha fact that resolution is sacriflcad for computational spaad. Thraa dl- 
manslonal models are bulky and time consuming but have much better 
resolution, however, when long term simulations are necessary, a one- 
dlmenslonal model Is recommended. 

The model described here can be modified to Include the single 
effects of the various quantities involved in the surface heat transfer 
phenomenon rather than using the equilibrium temperature concept* 
This is particularly recommended for the user who is interested in 
modeling the long term effects of one (for example, evaporation) of 
the quantities involved In the surface heat transfer processes.' 

Furthermore, jthe moctel can be easily adapted to handle connected 
multiple domain*. This recommendation is discussed in the text. 



SECTION 3 
PROGRAM DESCRIPTION AND PLOW CHART 

DESCRIPTION OF PROGRAM ALGORITHM 
Background 

A view of an idealized deep body of water is shown in Figure 1* 
This basin is divided Into eleven slices. The inner nine slices are of 
aqual thickness, 01, while the top and bottom slices are of thickness 
DZ/2. The thickness, DZ, is determined from the depth of the basin 
and the number of slices used. The temperature of each slice Is as 
shown In Flgurw 1 ; the horizontal lines correspond to the center of each 
slice. 

The condenser cooling water (CCW), If any, could be taken from 
any slice.. In Figure 1, the CCW is extracted from the center of Slice 2 
which Is at temperature T-. The discharge temperature, T Q , is the sum 
of T, and the increase in "temperature through the condenser. T Q is In- 
jected into a slice of equal temperature or treated as a surface outfall if 
T D is greater than the highest temperature of the basin. 

The basin also gains or loses heel from the surface as a result of 
changing climatic conditions which are required as input data. These 
could vary every time step, daily or monthly. 

Algorithm 

The problem is an initial value problem, so the values of dependent 
variables are assumed known Initially. The governing and associated 
equations are discussed In the next section. The governing equation is 
parabolic and mathematically represents a diffusion process with vertical 
convection. 

The values of the dependent variables at successive time steps are 
obtained by using a forward-time Dufort-Frankel scheme. 

The sequence in which calculations are performed is as follows: 
(Refer to Summary of Variables - next section.) 



The dependent variables, T, K,, W*, Ay, p, T~ and K», 
tialized. The area of each slice is calculated ana then trie 
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Figure 1. Idealized deep body of water 



Is calculated. .The heading of the beginning year Is prlntad. Tha 
valuas of tha variables, K-, W*, A v , p, T- and K-, are than calcu- 
latad. Tha temperatures 6T tha slices ara finally calculated* If tha 
tamparatura proflla li unstable, mixing of tha unstable portion of 
tha profile is undartakan. 

2. During the. next tima step, tha tamparaturas ara updated, and tha 
dependent variablas ara calculated again. 

3. Tha valuas of tha tamparatura T, addy diffusivity K z , numbar of 
days and surfaca haat transfer coefficient K- are prfhted every time 
step, every day or normally at the end of each month. At the end 
of the present year, the title of the new year is printed and compu- 
tions continue as listed above. These steps are shown in a flow chart, 
Figure 2. The results are stored on a magnetic tape and plotted 
when necessary. 

Description of Main and Subprograms 

The fortran calculation programs consist of a main program (NASA) 
and seven subroutines (YEARS, EQUIL1, STORE , CCW, SMOOTH, MIXIT 
and AREAS). 

1. MAIN; The main program handles the input data, calls the subrou- 
tines and does the temperature calculations. Two alternatives are 
given for handling the input data; these are either read through 
cards or in-data files or through a block-data arrangement given 

at the beginning of the main program. For users interested in the 
block-data package, the following caution is necessary: Whenever 
a data or set of data is changed, the main program must be recom- 
piled! 

2. YEARS: This subroutine prints the year heading. It is called at 
the beginning of a new year. 

3. EQUIL1 This subroutine reads the dew point temperature, wind 
speed and solar radiation. It then computes the surface heat trans<- 
fer coefficient and the equilibrium temperature. Depending on how 
the data has been averaged (e.g. days,, months or years); it is 
called as often as needed. 

4. STORE: This subroutine stores the calculated data on magnetic 
tape designated as Unit 8. The stored data could be read by the 
plotting subroutine called READER. This subroutine and other plot 
programs are described later. 

5. CCW: This subroutine supplies the condenser cooling water data. 
The data is also converted to the required units by this subroutine. 

6. SMOOTH: This subroutine finds the largest value of the eddy dlf- 
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Figure 2. Flow chart (calculation) 



fuslvlty and uses it to calculate tha variable time step. It also 
smoothens the calculated eddy dlffyslvity for unstable temperature 
gradients. It is called every time step. 

7. M1XIT; This subroutine looks for unstable temperature gradients 
and mixes or stabilizes the temperatures. It is also called every 
time rtep. 

8. AREAS: This subroutine handles the surface areas of each slice 
and converts the values to the required units. It is called only 
once at the beginning of the computations. 

9. INPUT: This is an in-data element containing all input data. 



SECTION <t 
DESCRIPTION OF PROGRAM SYMBOLS 



Introduction 



The programs have been written to calculate, as a function of depth, 
thermal diffusivity and temperature profiles over complete annual cycles. 
The equation integrated is 

A (Z)Tt<> c P T) 'W* ? A tzi K zlfr - fe«P c p A C« TV z ) +QA ' + \Z) H (Z) 

(1) 

The above equation requires two boundary conditions and one ititial con- 
dition. 

The initial condition is an input quantity supplied by the user and 
equals the homot hernia I temperature of the basin. The boundary condi- 
tions are: 



1. At the surface; 



K Z3Z|Z=h = K S (T E " T S ) ' (2) 



w.here Z = vertical coordinate measured from the deepest point 
Tg - equilibrium temperature 
T- = surface temperature 
Kg = surface heat exchange coefficient 

2. At the bottom; 

Perfect insulation is assumed, 

H-ii-o-* (3) 

Calculations of the temperature profiles are made by numerical integration 
of Equation (1). Calculations start with the homothermal conditions and 
a forward explicit scheme is used. 

Each time step, the surface temperature, T s = T12, is calculated 



and then the temperature of each slice is calculated. Solar radiation is 
absorbed at the surface slice and the unabsorbed portion is transmitted 
exponentially to the slices below. 

The empirical relations involved in this manual are summarized below. 
A full discussion is given in the final report, Lee et al. (1980). 

Description of Main Variables 

1. Density, p, fortran variable - ROW: 

P=A 1 +B 1 T+ C t T l (4,1 

where A, = density at 0°C 

1 = 1.02943 gm/cc 

B, = constant 

1 = -0.00002 

C- = constant 

1 = -0.0000048 

2. Eddy diffusivity, K z , fortran variable = XKZ 

K z = K ZQ (1 + ^R.J-1 (5) 

and 

R . *V g Z 2 3T m 



where R. = Richardson number 

cr! =0.1, an empirical constant, fortran variable - SIGMA 
g 1 = acceleration due to gravity, fortran variable - G 
W*= friction velocity, fortran variable - FRVEL 

= (T S /P) 

ay = A 2 + B 2 (T - 4) + C 2 (T - 4) 2 (6a) 

fortran variable for cu,, AV 

where A, = 0, volumetric coefficient of expansion at 4°C, fortran 
variable - A1 
B- = constant, fortran variable - A 2 

» 1.538 x lO"* 
C, = constant, fortran variable - A3 
* 3 -2.037 x 10 -7 

ayj can also be estimated by using Equation (4). 
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where K 7n = eddy dlffuslvlty under neutral condition (varies with 
iU time), fortran variable - XKZO 

K ZQ * A 3 + B 3 »ln(-||gt + C 3 ) (t l» in days) (6b) 

where A- * avaraga valua of K 7or fbrtran variable - R9 

Bg « half annual varlatloTTof K- , fbrtran varlabla - R10 
cf-3 phase angla, foreran variaTsfe - R8 

3. Haat source, H, fortran variable - F6 

H * n(1 - 0)A (Z j* o exp(-n(Z - h)) (7) 

where »0.5, fraction of the solar radiation absorbed at the surface 
n = 0.75, solar radiation absorption coefficient 
t = net solar radiation reaching the water surface (Input 
. ° variable), fortran variable - HSOL 
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SECTION 5 
PREPARATION OF INPUT DATA 



The Input data is storad In an in-data file - INPUT. Alternatively, 
it could be punched on cards. The Input data Is read in with an open 
format. The main variables read are; dewpolnt temperature, wind speed 
and solar radiation. In some cases where the dewpolnt temperature is 
not available, the relative humidity, air temperature and a pschometric 
chart are used to find the dewpolnt temperature. If this involves a lot 
of chart reading, subroutine EQUIL1 could be modified and the dewpolnt 
temperature calculated from a known oquation supplied by the user. 
If the latter case is used, then the input data base is enlarged to read 
air temperature, relative humidity, wind speed and solar radiation. A 
detailed input list of the constants is given in Appendix A. * 
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SECTION 6 
PLOTTING PROGRAMS AND EXECUTION ELEMENTS 



DESCRIPTION OF PROGRAMS 

The fortran plotting routine consists of one main program (PLOTTER) 
and one subroutine (READER). 

PLOTTER: This program calls the calcomp fortran subroutines (re- 
fer to a Calcomp plotting manual for details) and the subroutine (READER) 
which reads the calculated results from a magnetic tape designated as 
Unit 8. '(See Item A. 4.) A flow chart is shown in Figure 3. 

READER: Reads the calculated data stored on Unit 8 (magnetic 
tape) . 

Execution Elements 

Two execution elements are used, one for executing the calculated 
results and the other for executing the plots. 

DO-IT: This element compiles and prints the main program (NASA) 
and then prepares an entry point table, maps the necessary programs 
and subprograms, calls the In-data element containing the input data 
and finally, executes the calculations. This is done as follows for a 
UN I VAC 1100 computer at the University of Miami. 

Only one magnetic tape is necessary. 

1. ASC, AX FILE. 

The 'FILE* is assigned for the run. 

2. ASC, T 8., 16N, TAPENAME 

A magnetic tape file named '8.' is being assigned. The tape is 9* 
track, and the reel number Is TAPENAME'. The calculated results 
are stored on this tape. 

3. « PRT, S FIL2.NASA 

The main program is printed. 
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4. 9 PACK FILE. 

Th« 'FILE* is packed, 

5. PREP FILE. 

The entry point table is prepared. 

6. Q MAP, S 

7. IN FILE. NASA 

8. LIB FILE. 

9. END 

10. 9 XQT 

11. 9 ADO FILE. INPUT 

12. 9 FIN 

* 

PLOT-IT; Similar to DO-IT, bu« handles the plotting executions. 
For a UNIVAC 1100 computer the following cards are necessary. Two 
magnetic tapes are necessary. 

1. 9 ASG, AX FILE. 

2. 3 ASG, T 8., 16N, TAPENAME 

3. ASG, T 11., 16, PLOTTAPE 

A magnetic tape file named '11.' is being assigned. The tape is 7- 
track, and the reel number is 'PLOTTAPE'. The plots are stored on 
this tape. 

4. 9 PRT, S FILE. PLOTTER 
The plot program is printed. 

5. 9 PACK FILE. 

6. 9 PREP FILE. 

7. 9 MAP, S 

8. IN FILE. PLOTTER 

9. LIB FILE. 
10, END 



14 



11. 9 XQT 

12. « ADD FILE. INPUT 

13. e FIN 
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Figure 3. Flow chart (plots) 
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APPENDIX A 
EXAMPLE PROBLEM 

Th« model described In this manual was verified using monthly- 
averaged data supplied by Duke Power Company for Lake Keowee, South 
Carolina* Accordingly, the data discussed below apply to Lake Keowee. 

SITE DESCRIPTION 

Lake Keowee Is located 10 km west of Greenville, South Carolina. 
It is the source of cooling water for Oconee Nuclear Station (ONS), It 
was formed from 1968 through 1971 by damming the Little and Keowee 
rivers. A connecting canal (maximum depth 30. S m) [oins the two main 
arms of the lake. Flow out of the lake is through the Keowee Hydro 
Station. Lake Keowee also exchanges water with Lake Jocassee-pumped 
storage station. The three-unit ONS with a net capacity of 2500 Mwe 
started operating In July 1973. ONS operated on annual gross thermal 
capacity factors of 11, 28, 69 and 59% in the years 1973 through 1976, 
respectively i From 1977 to 1979 the factors varied from 65 to 75%. A 
map showing the geometry of the lake is given In Figure 1. 

PROBLEM STATEMENT 

Calculation of Parameters and Input Data 

1. The fortran variable DM (I, J) is a two-dimensional array containing 
the temperatures at the connecting channel between Lake Keowee and 
the Jocassee-pumped storage station. The data is averaged monthly. 
The units are In degrees Celciv- (°C). I Is the year counter and J 
is the month counter. The inp.»s for the first year are punched on 
the first card, the next year on the next card, and so on. Accord- 
ingly, each card contains twelve inputs in open format (real floating 
point numbers). 

2. The following fortran variables /constants are also read in with open 
format, five on one card. 

ZYEAR: starting year - 1971 (could be changed). 

DZ: thickness of an fnner slice (ft) - (maximum depth of lake)/(10.0) 

XKZL: lower limit of the eddy diffusivlty (ft 2 /day) - corresponds to 
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Figure 4. Lake Keowee 



the thermal diffusivity of solid water (15 ft 2 /day). 

H: maximum depth of lake, ft (150 ft). 

G: acceleration due to gravity (ft/sec 2 ) • 

PI: it = 3.1413926. 

A1: corresponds to A2 in Equation (6a); Al = °C~ . 

A2: corresponds to B2 In Equation (6a); A2 = 1.538 x 10" 5 °C. 

A3: corresponds to C2 in Equation (6a); A3 = -2.037 x 10 °C. 

A4: corresponds to A1 in Equation (4); A4 = 1.02943 gm/cc °C. 

A5: corresponds to B1 in Equation (4); A5 = 0.00002 gm/cc °C. 

A6: .corresponds to CI in Equation (4); A6 = -0.0000048 gm/cc °C 2 . 

(NOTE: The units for A4 through A6 are automatically converted to 
consistent units in the main program.) 

TO: homothermal temperature of lake (initial condition); TO = 7.8 °C. 

C i specific heat; C„ = 1.8 BTU/lb °C. 
P P 

SIGMA: see Equation (5); SIGMA = ^ - 0.1. 

**R6,R7,R8: the friction velocities (t /p) are calculated for the whole 
period and fitted into a sine curve: (friction velocity 
OMEGA) 

W* = R6 + R7 sin(-|~time + R8) 

where R6 = average value of W*, 0.1 ft /sec. 

R7 a average value of the half annual variations 

of W\ 0.025 ft /sec. 
R8 = phase angle, 2.61 radians 
TIME is in days, not specified. 

R8,R9,R10: correspond to C3, A3, and B 3 of Equation (6b) respec- 
tively; R9 = 800 ft 2 /day and R10 = 200 ft 2 /day. 

DATA1: or 1 (see below). 
3. The next set of inputs is the dewpoint temperatures, wind. speed and 

••Alternatively, friction velocity could be read in as monthly averages. 
If this alternative is followed, then DATA1 = 1, otherwise DATA1 = 0. 
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solar radiation. These can either b« punched on cards or stored in 
an In-data element. They are read every month. Each card contains 
three members. For example: for January-March 1971 (Lake Keowee) , 
the data are 

3.0, 6.69, 167.0 

0., 9.3, 264.4 

(i. 3, 9.28, 264.4 

The first number on each line (each card) is the dew point tempera- 
ture in °C. The second one is the wind speed in ft/sec. The third 
quantity is the solar radiation in BTU/ft 2 day. If DATA1 = 1, a fourth 
number must be included on each line (every card) . This fourth 
quantity is the computed friction velocity for each month. 

NOTE: The in-data element described above is called INPUT,. (See 
Fortran Source Program Listing, Appendix B.) 
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Sample Output and Sample Plots 
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■ • 

.3 .a .0 *a -a .o .a -a ,«j -a .c .it 

•nUNtH IS JAN. 1971 

3D ilia. * 1-IiO 7.dC 7.*C 7.80 7. SO - *-»w 7.80 „ 7.80 7. 60 l.*C 7.8C 

SC1.7J 801.70 401.70 801.70 SCl.lC 8Q1.JQ 801.70 831.70 «C1.70 8gl ,7n 8C1.7Q 801.70 4C1.70 

.3 .= .c .0 .Q .0 .0 .o .n .a .a .a 

•*QK|h IS fE9. 1971^ 

b r 7!i« l!»9 7.J9 7. iS 7.39 7.«9 7. IS 7.39 7.19 7. 19 7.39 7.i9 

JC1.SI 1C1-.S1 701. S7 701.57 7CI.S1 701.51 731.57 131.57 7C1 .57 lQl.5? 7CI.57 701.51 701.57 

,c .5 .c .a .c .o .3 .c .3 .w. .o .a 

10NIM IS r*»C« 1971 

SO *7*59 1- 7lfcC l.oC l.fc2 7.6to 7.77 7.95 6.1* « -I* »•'•* „ »;Sfc S.tS 

627. tC 570. £3 111.23 ;7L.fci 159.76 97.70 13.Ha lb. 81 tH.oi 221. ?Q 116.75 591.37 d27.bC 

.3 .3 .0 .0 .0 .0 .3 — - .0 • -0 .C .0 .0 
."OMH IS _»PHIL 1911 

12C "llll l *lC71 7.7i 7.60 7.52 t» . 11 8.SS 9.25 1C.52 11.22 11. VS Il.t" 

ccc.c: u.is ie.is ie.is ie.*s ia.15 ia.«£ is. is ia„is 13.91 12*. 19 m.si t-a.ci 

.a .3 .0 .0 .a .0 .0 .0 -3 -a .0 .0 

.10!. tn IS _«*» 1911 

15C *lll9 lB ll7* 1.43 1.S3 8.12 8.50 9.2 l lalsi 12.25 11. 65 *»;** »•?'* „ 

625. "0 IS. CO IS. 00 15.00 15. 00 15.0C 15.30 15.00 15.00 lS.00 19.22 3U.56 625.90 

.a .c .o .o .o >o >c «o , .o .o .0 .0 

- JlOhlb IS JUNE 1971 l 

1*C *7*92 K92 7.99 8.15 8.17 9.11 10.31 1Z.50 1*. 2* 22.21 « 2 '»» „ 23.16 

698.61 15. CO IS. 00 IS. CO 15. CC * lS.fiG 15. CG 15.00 15.00 IS. 00 17.12 11C.3S 696.61 

.0 .0 .0 .C .0 .0 .0 .0 .9 .0 .0 .3 

y HOMIh IS JUL¥ 1921 ...... ...— .. 

* 210 21.15 2C1.11 . „ . 

210 • .!«• 8.M 8.2* 8.51 9.0J IC.G1 11.77 11. 11 19.29 23.83 '*•»' .» 2 *;12 ,. 

798.29 IS.CO 15.00 15. CO 15. OC 15.00 15. CO IS. 00 15.0a lS.CO 32.60„ 622.68 198.29 

.q .a .o .o .c .o .o . .a .a .0 .0 «o 
•0N1H is »ur. . 1971 

210 2 %li» U 8^8 8. 61 9.Ci 9. 71 11.15 13. iC 16.93 22.1'» 21.21 2*. It 2l.S« 

898.12 15. CO IS. 00 li.rC 15. OC 15. &0 15. CG IS. 03 15.r0 15.00 S3. 53 771. H 89*. i2 

.C .3 .0 .C .C .C .0 -3 .0 .0 .0 .0 
-Ot.fh IS SEP I. 1971 

273 'e'98 BVIB 9.2C 9.7« 13.75 12.12 l«i.n7 I£.-«i 22-6T 22.98 23. Gl 23.12 

972. 2* 15. CO 15.00 15. ZL 15- CC 15. -Q 15. CC 15.30 15.30 ;e . «5 12S.C9 926.51 912.39 

.c ,c .o .o .o .e .n .n .Q .c .o .3 

-OMIn IS OCT. 1971 

32C B.21 167.98 

3Cr 9.63 9.63 «.*; It. 56 1 1 . »«► 13.4.7 16. is |<>.i5 19.71 IV. 76 19.77 19.17 

999.99 15. CO 15.00 15. CO 15. OC 15.00 IS.BC 1 S-3D 20 . 8 1 J-t.77 8-71.99 996.95 999.99 

.3 .o .c .0 .0 .o .a .o .c .0 .o ,o 

-OhlH IS NOV. 1971 

3 J" 8.73 16 8.76 

33C iC*2 10.12 10.71 11. ie , 12. li 11.21 1«.23 11.23 H.23 11.23 11.23 11.25 

971.11 15.C0 IS.OQ 15.00 15.00 15.00 23.16 83C.29 971. 11 97i.il 91i.U 91l.ll 971. II 

.0 .0 .0 .0 .0 .0 .0 .3 .0 .0 .0 „0 

nQHIH IS DEC. 1911 

36C 3.59 132.35 ,. ., „ 

36C 11.96 11.17 11.58 11.95 11.9* 11.96 11.96 11.96 11.96 11.96 11.46 » l **f „ 

901.10 171.92 92.91 3C.18 18.01 114.19 931. 1C 901.10 9CI.1C 9Q1.-1C 931.10 901. 13 901.10 

Figure 5. Sample output - Lake Keowee, 1971 



TEMPERATURE PROFILES FOR LRKE KEQWEE 
< DEPTH IS HERSUR2D fRQrt THE DEEPEST TO INT OF THE LRKE J 
(STATIONS SQQ-506) 
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Fiaure 6. Sample plots - measured average temperature profiles (Stations 
500-506) vs predicted temperature profiles, Lake Keowee, 1971 
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STRRTIFICRT1QN CYCLE FOR LRKE KEGKtE 1971-i9T9 



o 
o 



to 



Solid Lines (No Discharge) 

Broken Lines (Discharge - Mid-layer Temperatures) 



^equilibrium Temperatures (*} 

Surface Temperatures^) 

n 



Mid-layer Temperatures (a J 




^'.00 



36-00 7Z.GQ 108.00 144.00 

TIME IN DfiYS kIO 1 



lSlfi 1S76 .11 1S11 1978 1979 

180-00 zAj-QO ZSZ.OO Zflfi.OO 3Z4.Q0 



Figure 7. Sample plot 
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FORTRAN PROGRAM LISTING 
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NASA SYM CREATED ON 12 AuG 80 AT 1«*:17;D;> 

1 C ONE DIMENSIONAL MODEL FOR THE -SEASONAL THERMOCLINE 

2 C 

3 C 

H DIMENSION T (20 1 ,AV 1201 ,CB ( 20) . 2( 20J , A I 2D) , XKZI 20 I ,ROU 1 20 J ,TN 1 201 

5 DIMENSION DM (20) , T2(201 ,XTDO( 1 0, 36Q1 

6 DIMENSION LELT EM( 1 21 . QP( 12 1 

7 CHARACTER*6 MONTHS (12) 
ft C 

9 C 

10 DATA(H0NTHS(J1 ,J=1 ,121 /'JAN.* , •FEB.", •MARCH', • APRIL 1 , 'MAY • , 'JUNE 1 , 

11 C'JULY" .'AUG. •.•SEPT. •.•OCT. • , • NOV. • , 'OLC. • / 

13 L 

1" C IF YOJ NEED 10 STOKE RESJL1S ON MAGNETIC TAPt READ JRJN=1 

15 C OTHERWISE JRUN-2. 

16 C 

17 READ ] ,JRUW 

IS REAO 1 ,IYEAR,0Z,XKZL,H,G 

19 READ 1 ,PI, A1,A2,A3,A<| 

£ *Ii READ 1 e AS .Afc.TO.CP, SIGMA 

21 READ 1 ,R6 t Rl,R8 I R9 l R10 

22 1 FORMAT (1 

23 HMI=0 

2*» zm-o. 

25 JIM=1 

26 TDD-O. 

27 DVE=C. 

28 CALL AREAS (A) 

29 j.=i 
3D JUri 

31 JJ-0 

32 NDAYS=Q 

33 NDAYS1=0 

34 TIME=0. 

35 T1ME1-U. 

36 TIME?=0.. 

37 TlME3rQ. 

38 T1ME4-0. 

39 TE-TO 

40 DO 20 1=1 ,12 

41 TIII=T0 

42 T2U1-T0 

43 20 CON I IN JE 

44 DO 22 1-2,11 

45 2 fl l=DZ/2.*(l-2 1*D2 



46 22 CONTINUE 

47 Z(12l-H 

48 DT=(0.4*0Z**2l/100U.U 

49 GP2=57l.G73e3*(6a.**2»*21 . 

50 CALL YtARS (SELTEM, uOPP.I YE AR1 

51 CALL CCUIQP.DELTEM.IYEAR.DTl 

52 U-Q 

53 OMEGA =2.*^ 1/360. 
51 T(l2l=T0 

55 T12--T0 

56 JTOT-1 

57 MJ=1 

56 ROW (12 1 = ( A1 ^5*1(1 21 *A6*lT( 121 )** 21*62.1 

59 R0kCP=R0U(l2l*CP 

60 CALL EttUlLl (IN a T£t 4K.TDE.tr ■*IND»HSGU 

61 33 IFIHJ.EQ.l 1DELTm2=oM(1 l~T(7l 

62 FRVEL=(R6 *R7*S IMD HEGA*TI HE*R8 1 1**2 

63 XKZ0=(R9*Rl0*SINf0KEGA*TIHE*R8iI 

61 AVC1 1=A1*A2* tTIll -*l. 1*A3*«T«1 l-*i.l**2 

65 XKZ U JzXKZO* (1»S1GHA*A\|!( 1)*G*IIH-2I 11 1**21* 

66 1(3.*T(1 1*T(31-1.*T(211/I2.*DZ*FRVEL11**CN-11 

67 00 90 1=2,11 

68 AV(11=AI*A2* IT (II -1.1* A3* (Til 1-1.1**2 

k, 69 X*2 II 1=XKZ0*(1*SI6KA*AW(I1*G*I(H-Z(11 1**21* 

«» 70 1 (1 (1*1 1-1 41-11 1/(D2*FRVEL1 l**(N-ll 

71 ROW (I l=(A1«A5vT lI)«A6*f Till 1**21*62.4 

72 90 CONTINUE 

73 R0WI12 1 = ( Aq*A5*T I 1 21*A6* ( T ( 121 1**21*62.1 
71 AV (121=A1*A2* (1(121-1 .1*A3*(T( 121-4.1**2 

75 XKZ1121=XKZ0*(1*SI6MA*AV(121*G*S (H-Zl 1 21 1 **21 * 

76 I I3.*T(11 1«T (91-4.*1 HOI 1 /(I.5*0Z*FR\/ELl 1**|N-T 1 

77 R0wCP=R0W(121*CP , 

78 CALL iM00ThlXKZ,XKZU l XKZL,hDAYSl,TM2, 112,T,DI1,0Z1 

79 902 DO 989 1=1 ,12 

bO IF (XKZ (I 1.LT.XKZL1XKZ(II = XKZL 

111 IF (XKZlI l.LT.XKZ01XKZIIl=XKZ0 

62 989 CONTINUE 

B3 DO 91 1=2 ,11 

81 F1=DT/ (ROW (11*CP*A(111 

85 F2=( (R0U(Il+R0w(I«ll}/2.*(A(Il*A(l+lll/2. 

86 1*IXK?II 1*XKZU*111 /2.* (1(1*11 -1(11 l-( (ROwlIl 
fc7 2* ROW (1-1 1 1*(A (Il^A (1*1 11/1 .*(XKZ(I1 

68 3*XKZ (1-1 J 1/2.*(T( Il-T( 1-11 111 /(DZ**21 

69 IF (lYEAR.Lt.l973lDt_LTH2=0.G 
^0 IF(IYEAR.LE.1973ICP2=0.0 

91 777 F3=RGw (I 1*L£LTEH( J»l*CP*liP(JUl 

92 F31=R0*(I 1*0ELTM2*LP*QP2/A(I1 

93 F11=(R0U(I1*CP*GP2/(1.5*DZ)1*DELTH2*(T(I*11~TC1~1 11 

*1 Fl = (RO«r(Il*CP*QP( J*J/(1.5*0Z1 1 *0ELTEH(JU1.*(T(I*1 1-T (I~l 1 1 

95 IF1T(I*11.LE.T(I-1 1 lF1 = l!R0ii (1 1 *CP*tf> ( Jkl/ ( 1 .5 *GZS 1*DELTEH(JU|' 



96 IF IT il *1 l.LE.TlI-l 1 IF* 1=«R0UI 1 l*CP*QP2/f 1 . b*D2l 1*D£LTM2 

97 F5-TII1 

98 F6=0.5*lEXPI-0.75*tH-2tII) ))*<HS6L) 

99 F7=-0.75*Aill 
100. Ffc--F6*F7 

101 TU-TI8 >«DELTE.M(JU> 

102 IF (I.LT.8 IXAK=U. 

103 1FfI.GE.BlXAK=l. 

104 IFITIII.GT.TOlXKrG. 

105 JFCTIII.L£.T01Xn=l. 

106 TD2=0ELTH2«T<5I 

107 IFC1 I] ).GT .TD2)XM1=0. 

108 IFITII I.LE.TD21XM1-1. 

109 IFU.LE.S)XTK=1. 
1IQ IFlI.GT.5»XTK=l. 

111 TN<II = (F24F3«XAK*XM*F<i+xril+Fm*F31*XTh*F8)*Fl*F5 

112 91 CONlINJE 

113 TNU)=TI21 

111 TH = fTN't'12)«T'DEWI/2.0 

115 Fk=9.2*0.16»CUlHD**2l 

116 BETA=0.35*0i015*TH*0.0012*lTH*»2l 

117 XK=Ct. S+0.Q5*TN(121+eETA*Fri*0.47*Fhl*<l. 232*19. /5. J 

118 TEzl[|£W*HSOL/XK 

119 C0NS1-I1 .S*XK*D2) /(R0bCP*XK2< 121 J 
w 120 TEH =11411 11 

o 121 T£10=TN(10> 

122 SHEAT=(ROUCP*DELTEKC JU»»UPIJM) )/(A(12l*XK) 

123 TFflO.GT.lNtl2))Gu TO 11 
121 GO TO 15 

125 11 TKI12 )=<1.*TN(11) -THI 10) +C0NS 1*1E*5HEAT*C0NS1 1/C3. ♦COhSl 1 

126 CO TO 16 

127 15 TUI12 > = (<». *TN(ll>-?N(ia)+C0NSl*TEl/l3.»C0NSll 

128 16 TS=TNI12) 

129 CALL HIXTT UN ,A1 

130 TIHE-T1ME»0T 
121 TIME2=TIME2*DT 

1 32 TIME3=TIHE3«0T 

133 TIME1=TIME1*DT 

131 TlME5rTlMES*DI 

135 DO 929 Trl ,12 

136 N T2IHrlN*(II 

137 929 CONTlNUC 

138 T12=TI12» 

139 600 TK12-=TNI12 1 
1*40 601 00 92 1-1 ,12 
111 TlIlrTNII) 
142 92 CONTINUE 
113 J=J*1 

1H TIME1-TIME1*0T 

115 IF If.OAYS.GE. 360 1TIME3 = 1I HE 3-3fc0.0 



1«6 1MNDAYS.GE.3601T IME2-TI ML2-360 . O 

ikl IF INDAYS.CE.360)TIHL = TIM£-36 0.0 

1H8 IFCNDAYS.GE.360)T 1KL4-T1 Mt<*-360. & 

1*9 " IF <NUAYS.GE.36C)TlMt5=TlME5-360.D 

150 IF<NDAYS.GE.360)JJ=0 

151 IF«NDAYS.bE.360iJ.rj 

152 IKIYEAR. 67. 19791 60 TO 99 

153 IFrN0AYS.GE.360) IYEAR=IyEAR*1 
ISil IF (NDAYS.6E. J601CALL CCW « uP.DEL TEH, I YEAR, DT J 

155 IF(NDAYS.GE.3601CALL YEARS! SELTEM ,QCPP,I YtAR I 

156 IF Cf.DAYS.GE.360J jr0T = JT0T«l 

157 IFINDAYS.GE. 360) JIM=JlM*l 

158 IF(TlME1.6E.1.0)GC TO 501 

159 GO 10 &02 
ILO 501 NHI=HHI«1 

161 XTOD|JlH t HMIl=TD 

162 TIME«<-TIMEM-1. 

163 502 CONTINUE 

lt«« IF|NDAYS.GE.36Q)NDAYS=G 

165 DO 66 1:2,10 

1G6 COII 1-IT (Nll-TlIJ 1/15 . 

167 66 CONTINUE 

168 Cbll )=CT*2>-T 11)1/7.5 

169 CBI11 ) = (T(12I-If UJI/7.5 

170 IF CTIHE1.GE.30.) 6a TO 98 

171 TDD=TDD«TD 

172 OVt=DVE«l.. 

173 60 TO 3i 
171 98 NDAYS=TIME2 

175 TDL^TOU/DVE 

176 PRINT 988 ,ICP UrfJl f JkJ=I f 12) 

177 988 FORMAT UX,12F10.1i 

178 TlME«t=0. 

179 MHJ=G 

180 JJ-JJ-*! 
lfcl JU=J««1 
lte.2 N0AYS1=TIHE3 
lfcl HJrMJ^l 

181 0ELTK2=DH fKJ)-TC5) 
185 IF IMJ. uE.12)HJ=i 
le6 313 CONTINUE 
187 DO 700 1=1 ,12 
lfc8 710 TII1=T'N|I| 
169 IFURUN.-EO.ZIGO TO 111 
190 C - - - - 

1*1 CCI _ 

J »? PK^0.TE f MDAY"S,fN12 f Ti2,fl,F2.F3;F3l',F^lVF5TF6VFKF6TfD,TD2 

i?> cndaysi.tihe.iimez, hmes.Iyear.hj.xkJtog.j) 



FCJRUN.E0.2) GO TO 111 

:alL 5T0REIT,AV,CB,2 t A,XK2 f R0U f TN, DM . T i f MONTHS ,T2 .GP , 
:P 1 SIGMA t Ri f Rq.R5,k6 I R7,R8,R9 t Rl0 fW P2 f FREVtL l R0«CP i Df , 

K/:0 t TE f MDAYS,fN12,Tl2,fl,F2.F3 f F3l,F*Jl,F5 1 F6 t F7,F£,TD, 



191 


111 


195 




196 




197 


920 


198 




199 


10] 


2 GO 




2G1 




202 




203 




2 04 


18 


205 




206 


12 


207 


9 


208 


7 


209 




210 




211 




212 




213 




21«l 


99 


215 


921 


216 




217 




218 





CONTINUE 

CALL EQUIL MTN.TE.XK.TDEW .WIND ,H SOLI 

PRINT 920 .MONTHS < J J), f YEAR 

FORMAT <2X , 'MONTH IS«,2X, A6 , 2X , 111 

PRINT 101 ,N0AYS,TE .XK 

FORMAT UX,I6,2F9. 21 

UR1TEI6.91 NOAYS.C Till, 1=1, 121 

URITE(6,71 XKZO.CXKZI 11,1 = 1,121 

IFt(IYLAR.ES*l973.AN0.NDAY^.GE .21G1.0R. < I Yt.AR. GT. 1973 i J 
CURITEC6 .18 IT0Q,0EL1EH(JU-11 

FORMAT <1X ,'IHt AVERAGE MONTHLY OISCH. 1EMP. = • t F5.2,5X, 
C^ELTA-T = *,F5.21 

FORMAT C1X,11F10.21 

FOHHAT (IX .I6.12F9.2) 

FORMAT UX.13F9.21 

TlMEl=TlMEl-30.0 

TOOrO. 

DVE=0. 

IF IIYEAR.GT.197916ti TO 99 
60 TO 33 

PRINT 921 ,J 

FORMAT (2X. 'TOTAL NJMBER OF COMPJTATIONS = ',115^ X 12»J 

END FILE 8 

STOP 



w 218 END 



APEAS SYM CREATEC ON 12 AJG SO M 11;05:27 

1 C THIS SUBROUTINE CONTAINS \E AREAS OP 

2 C A DOMXN (LAKE KEO-EE), AT TWELVE 

3 C HOklZONTAL CROSS-SECTIONS. 
<♦ C 

5 , C 

6 SUBROUTINE AREASUI 

7 OIMENSION A(2Q> 

8 AC0NS=10.**8 

9 A (1)=0.032S*ACCNS 

10 A (2)=Q.055*AC0NS 

11 A m=Q.20C*AC0NS 

12 A («*)r0.55Q*ACONS 

13 A (5)=1.125*AC0NS 

14 A (b)=l.a*ACONS 

15 A <7I=2.575»AC0NS 

16 A (o l=3.5S*AC0NS 

17 A m=«*.7l>AC0NS 
13 A (10 )=S«82S*AC0NS 
19 A (11 )=7.2S*AC0NS 
ZO A(12 )=8.0Qd*AC0.\S 
21 RETJRN 

2 E.ND 



* m 

6, 



)b POOR QUATvm, 



QUAltTY 
33 



•CCW SYM CREATED ON 12 AUG 50 AT 13:00:09 

1 C THIS SUBROUTINE CONTAINS THE CONDENSER 

? C COOLING WATER. ASSUMES THAT COMPUTATIONS 

3 C START IN 1971. 

4 C 

5 C 

6 C 

7 SUBROUTINE CCW <QP, UELTEM, IYEA*,0T> 
3 DIMENSION (*Pfl2> ,0£LTEMf L2> 

9 IF (IYEAR.GT. 1979)00 TO 11 

10 It£A=IYEAR-1970 

11 AC0ST=10.0 

12 GO T0(i,l,3,4,S f 6,7»a f 91,lY£A 

13 1 DO 10 1=1 ,12 
Id QP(Il=U.O 

15 IP DELTEMU)=C.O 

16 GO TO 11 

17 3 DO 12 1=1, o 

ia op(i)=o.a ' 

19 12 OELTEMtIl=L.Q 

20 OP (7 1=1890 »2»ACQST 

21 QP18 1=1910. 3*AC0ST 
22- GPC9 1=2170. 7*AC0ST 

23 QP (1C1=2232.5*AC0ST 

24 QP(11)=217C.7*ACQST 

25 OP 112) =3284.6*AC0ST 

26 DELTEM (71 = 5.3 

27 OELTEM (3 1=4.6 
25 OELTEM (91=5.3 

29 DELTEM (10 1=7.3 

30 DELTEM (11 1=7.7 

31 OELTEM (12 1=4.1 

32 GO TO 11 

33 4 QP (1 )=3069.3*AC0ST 

34 QP (2 1=3069. 4*AC0ST 

35 QP (3 1=2976. 9*AC05T 

36 QP (4 1=2307. 3*AC0ST 

37 QP<5 1=2164. 6*AC0ST 

38 QP (6 )=4171.8*AC0ST 

39 0P(71=S334.6*AC0$T 

40 QP (8 1=4727. 1*AC0ST 

41 QP 19 1=5961. 4*AC0ST 

42 OP (10 1=4953. 4*AC0ST 

43 QP (11 1=4202. i*ACOST- 

44 QP (12 1=5225. 6*AC0ST 

45 OELTEM (1 1=U. 2 

46 DELTEM (2 1=7.4 

47 DELTEM 131=8.4 
43 OELTEM (4 1=8.0 

49 OELTEM (5 1 = 2.7 

50 OELTEM (61 = 6.0 

51 OELTEM (71=5.0 

52 DELTEM (81 = 4.8 

53 OELTEM (91=5.8 

54 OELTEM 110 ) = 3. 5 

55 OELTEM (11 1=7.9 

56 OELTEM (12 1=5.9 

57 GO TO 11 

58 5 QP (1 1=4612. 4*AC0ST 

59 QP (2 1 = 3694. 9*AC0ST 

60 QP (3 1=5456. 8*AC0ST 

61 QP (4 1=5570. 3*AC0ST 

62 QP (5 1=6494. 3*A COST 

63 QP (6 1=6574. 2*AC0ST 
b4 QP (7 1=7104. 2*AC0ST 
fe5 QP (8 1=7510. 1*AC0ST 

66 QP (9 1=7201. 6*AC0ST 

67 QP (101=6993. 4*AC0ST 



34 



T~;y?.#:?23?; : -r' rt - 



63 

69 
70 
71 
72 
73 
74 
75 
76 
77 
73 
79 
60 
8 1 
82 
83 
«4 
65 
36 
67 

aa 

89 
90 
91 
92 
93 
94 
95 
96 
97 

9a 

99 
ICO 
1C1 
1C2 
1C3 
1C4 
IDS 
1C6 
1C7 
108 
1C9 
110 
111 
112 
113 
114 
115 
116 
117 
113 
119 
120 
121 
122 
123 
124 
125 
126 
127 
123 
129 
130 
131 
122 
133 
124 
135 
136 
137 
1 3fl 
139 



OPll 
QPC1 
DELT 
DELT 
OELT 
OELT 
OELT 
OELT 
OELT 
OELT 
OELT 
DELT 
DELT 
DELT 
GO T 
OP (1 
OP 12 
QP(3 
OP (4 
CP (5 
OP (6 
OP (7 
QP (8 
QP(9 
OP U 
OPll 
OP 11 
DELT 
OELT 
DELT 
DELT 
DELT 
OELT 
DELT 
OELT 
OELT 
OELT 
OELT 
DELT 
GO T 
QP (1 
OP 12 
QP(3 
QP (4 
QP (5 
OP (6 
QP17 
QP {8 
OP (9 
QP (1 
QP 11 
QP 11 
DELT 
DELT 
DELT 
DELT 
OELT 
DELT 
DELT 
DELT 
OELT 
OELT 
DELT 
OELT 
GO T 
OP (1 
0PJ2 
QP (3 
OP (4 
OP (5 
QP16 
QP 17 



11=7146 
2) =68 5 
EM ( 1 ) = 
EM (2) = 
EM (3 1 = 
EM (4 ) = 
EM (5) = 
EM (6) = 
EM(7) = 
EM tfl) = 
EM(9) = 

EM (10) 
EM (11 I 
EM (12) 

11 
1=606 9 
1=4440 
1=4874 
>=4272 
1=3970 
1=5197 
)=5830 
)=724S 
1=678 5 
0)=563 
1 » =5 8 
2) =491 
EMU) = 
EM (21 = 
EM (3) = 
EM (4 ) = 
EM (5) = 
EM (6 ) = 
EM (7)= 
EM (8 ) = 
EM 19) = 
EM ( 1 ) 
EM (11 1 
EM (12 ) 
U 
1=5045 
)=4985 
1=5113 
)=6013 
1=6 302 
)=439S 
1=5033 
> =5 70 8 
1=6964 
1=6 75 
1 )=469 
2)-S8S 
EM (1 » = 
EM «2 ) = 
EM < 3 > = 
EM (4) = 
(5) = 
(6 1 = 
<7> = 
EM (8 1 = 
EM ( 9 ) = 
EM (10) 
EM (11 ) 
EM (12 ) 
11 
1=6176 
1=6444 
1=5195 
)=4811 
1=4984 
>=56 5 9 
1=70 53 



7.1*AC0ST 

C.9*AC0ST 

6.2 

4.o 

£.2 

c.3 

6.8 

b.6 

8.3 

7.8 

7.4 

=7.7 

=o • 5 

=9.4 

• 3*A 

.2*A 
.3*A 
.1*A 
.7*A 
.6*A 
.C*A 
,3*A 

• 4*A 
7.8* 
9.2* 
4.8* 
10.6 
7.3 
7.1 
5.1 
5.8 
9.3 
7.4 
o .*5 

b.o 

=7.8 
=6.7 
=6.4 



COST 
COST 
COST 
COST 
COST 
COST 
COST 
COST 
COST 
ACOS 
ACQS 
ACOS 



ST 

ST 

ST 

ST 

ST 

ST 

ST 

ST 

ST 

OST 

OST 

OST 



EM 
EM 
EM 



•8*AC0 

.2*AC0 

.5* A CO 

.6*AC0 

•4*AC0 

.3*AC0 

.6*AC0 

.9*AC0 

.0*ACO 

4.7*AC 

7.6*AC 

4.6*AC 

12.5 

11.4 

10.4 

11.4 

9.4 

6.4 

7.4 

5*0 

b.O 

= 3.8 

=6.2 

=7.9 



.7*AC0ST 
.6*AC0ST 
,7*AC0ST 
.8*AC0ST 
.2*AC0ST 
.9*AC0ST 
.6* A CO ST 
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HO 
141 
142 

HI 

144 

1*5 
146 
147 
148 
149 
150 
1S1 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
•164 
165 
166 
167 

ua 

169 

170 
171 

172 
173 

174 
175 

176 
177 
178 
179 
180 
181 
182 
163 



11 



0P(8 »=79l 
QP<9 >=655 
QP (10) =74 
QP(11)=60 
QP<12)=65 
DELTEM (1 > 
DELTEM (2> 
0ELTEM13) 
OELTEM (4 ) 
OELTEM (S > 
OELTEM (6 ) 
OELTEM (7) 
OELTEM (fl ) 
OELTEM (9) 
OELTEM (10 
OELTEM (11 
OELTEM (12 
GO TO 11 
QP (1 )=720 
OP (2 )=731 
QP (3 ):741 
QP(4)=727 
QP (5 )=418 
QP (6 >=S38 
QP(7 )=473 
QP(8 1=473 
QP(9 )=473 
OP (10)=47 
QPU1 )=4 7 
QP (12)=47 
DELTEM (1 ) 
OELTEM (2 I 
DELTEM<3 ) 
OELTEM (4) 
DELTEM 15) 
OELTEM (6) 
OELTEM (7) 
OELTEM (8) 
OELTEM (9) 
OELTEMUO 
OELTEMdl 
DELTEM (12 
RETURN 
E NO • 



4.9*AC0ST 

7.3»AC0ST 

Q7.4*AC0S 

65.1*AC0S 

GJ.5+AC0S 

= 9,0 

= 11.0 

= 13.2 

= 9.7 

= 1 C . 1 

= 6.1 

= 7.9 

= 7.5 

= 7.6 

)=6.2 

)=S.4 



7.7«A 
9. 9* A 
9.5 .* A 
5. 3w A 
9 . I * A 
i.2*A 
3.3*A 
3.3*A 
3.3*A 
33.3* 
33.3* 
33.3* 
= 10.3 
= 10.4 
= 9.6 
= 9.9 
:c.2 
= 7.1 
= 5.0 
= 5.0 
= 5.0 
1=5.0 
» =5.0 
) =5 . 



COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

ACOST 

ACOST 

ACOST 
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E0JIL1 SYM CREATED ON 11 JuN 60 AT 11:00:00 

1 SUaROOTINE E0UIL11 TN.TE , XK * TOE W,X TN, XTE, XXK.Ml ND, HSOL 1 

2 DIMENSION TK(20» .XTNI20I 

3 2 RE AD 15 ,1 ITLEW ,UlND ,HSOL 
1 1 FORMAT!) 

5 UIN0=bIND*U»4b 

6 HS0L=HS0L*3.6855 

7 TM=CTNU2 )*TDEU|/2.Q 

8 FU=9.2*0.«»6*JUlN0**2) 

9 BETA-G.35*0.015*TM*0.0012*(TM**2l 

10 XK=4.5*0.05*TN(121*BETA*FU*0.«»7** k 

11 XK=XK*<|. 232*19. /5. 1 

12 TE-TDEU+HSCL/XK 

13 XTM=IXTNI12J*T0EU) /2.0 
11 XF.< = 9.2*0.<*6*<WIND**2» 

15 XBETA=C.35*0.015*XTM*0 .Q012*C X TM**2) 

16 XXK=4.5*0.05*XTN(12l*XBETA*XFU*G.«*7*XFW 

17 XXK~XXK*<».232* I9./S.* 

18 XTE=TDEW»HSOL/XXK 

19 R£TURN 

20 ENO 



1 

2 
3 

1 
5 
6 

7 

a 

9 

10 
11 
12 
13 
14 
15 
16 
17 
13 
19 
20 
21 
22 
23 
21 
25 
26 
27 
23 
29 
30 
31 

33 
31 
35 
36 
37 
3S 
39 

«40 
11 

12 
13 
11 
15 
16 
17 
13 
19 
50 
51 
52 
53 
51 
55 
56 
57 
53 
59 
cO 
61 
62 
63 
61 
65 
66 
67 
63 
69 



INP 


JT S 


3.0 


,6.6 


C, 


9.3, 


6.3 


,9.2 


7.5 


.6.7 
2,7. 


17. 


IS. 


3,5. 


20. 


,6.1 


19. 


41t5 


13. 


33,5 


13. 


ae,7 


2 . 8 8 ,7. 


5.5 


, i,3 


1.67,6, 


-2. 


22,9 


1.11 ,9. 


b.67,8. 


11. 


11,7 


13. 


13,7 


13. 


77,6 


22. 


2 2,6 


18. 


3,5. 


a. 


5,7. 


5.9 


', 7.1 


1.1 


6.6, 


1., 


7.22 


-1. 


, 7.3 


la. 


r 7.1 


7.7 


,8.1 
3,6. 


11. 


20. 


25,3 


22. 


2,5. 


21. 


7,5. 


20. 


3,6. 


13. 


5,7. 


7.2 


,8.1 


3.2 


,5.6 


S.2 


,5.3 


c, 


5.3, 


6.3 


,7.7 


10. 


7,3. 


17. 


2 ,6. 


i?. 


8 ,6. 


21. 


,5.2 


21. 


,5.3 


17. 


5 ,6. 


10. 


2,5. 


6.0 


,7.2 



3.8,6.9 
3.0,6.3 
3.5,7.6 
2.2,9.6 
7.2, 7.6 
17.5,1. 
19.0,5. 
21.3,5. 
21.0,5. 
16.2,7. 
12.1,7. 
7.9, 
i y f 



6.9 

7.2 



-1.0 ,7. 
2.2 t 8.5 
3.9, 

11.2 

11.0 

13.3 

19.8 ,5. 

13.0,6. 

15.1,7. 



7.9 
,7. 



YM CREATED 


9,167.0 
261.1 


3,261.1 


2,157. 


5,130.5 


65 ,178. 


8,109. 


.75 ,123.2 


.77,329. 


.02,261.3 


53J217.7 


,117.7 


69 ,178. 


.26,257.6 


22,352.5 


72,118. 


.53 ,133.6 


.95 ,561.3 


.61 ,193.8 


.07,153.5 


17,336.3 


17,298.1 


3,220.9 


118. 


,162.7 


,279.5 


,318.5 


1,119.3 
83,119.5 


.01 ,5G7.7 


32,196.9 


1,391.6 


303,333.1 


1,311 .7 


1,217.6 
,151. 


,191.1 


226.9 


,326.1 
73,397.7 


3,13b. 


98 ,559.3 


,159.5 


7,180. 


71 ,339.2 


7,302.5 


,231.1 


,181.9 


93,191.1 


11,226.9 


,326.1 


,397.7 


8,136. 


32,559.3 


10,159.5 


1,180.3 


3,339*3 


7,302.5 


,231. 1 


,131.9 


1,209.3 


,310.9 


,333.6 


6 ,1'J6 .9 


3,118.1 


1,160.2 


9,188.3 


65 ,13C1 


13,315.1 



ON 12 AJG 80 AT 13:01:13 
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73 
71 
72 
73 
74 
75 
76 
77 
78 

79 

80 
81 




20.2,5.8,551.8 
20. 7, 5. 4, 423. 9 
IS. 7,5. 3,350.7 

82 9.2,7.2,256.6 

83 7.0,7.5,196.2 
64 0. 4, 7.2,178.2 
85 -2.9,7.9,227. 
8 6 -5.0 ,6.8,3 C8. 
bl 1.2,7.6,408. 

88 9.6,7.6,429. 

89 14. ,6.7,513. 

90 19.4,4.7,593. 

91 20.8,5.7,568. 
20.8,5.1,461. 
IS. 5 ,5.7i3b5. 
9.3,6.6 ,369. 
9.0,5.8 ,232. 
0.4, 7.3,191. 



92 
>3 

94 
"95 

96 _....._, 

97 -3.3 3,8.6,208. 

98 CO, 7.2,251. 
b.dr 7.9,373. 
S.2, 7.6 ,479. 
1 4 . f 6 . 7 , S 1 3 . 
i c *i a -r c a a 



99 

100 

1 U 1 !<•, t ci t f f 31ii 

1C2 19.4,4.7,598. 

103 20.8,5.7,568. 

104 20.3,5.1,461. 
1C5 15.5,5.7,387. 
106 9.3,6.6,36 9. 
IC7 " ~ ' " "-" 
103 
1Q9 



7. J , O .O , JO V • 

9.0 , 5.8 ,232. 
0.4 , 7.3 ,191. 
0.4, 7.3 ,191. 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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I 

2 

3 
<♦ 

s 

6 

7 
9 
9 

10 
11 
12 
II 
i«» 
15 
16 
17 
18 
19 
20 
21 
22 



MIXIT 

C 

C 

C 

C 

C 

C 



ICO 



SYH CREATEC ON 12 AUG 80 AT 13:26:57 



1 

10 



300 



THIS SUBROUTINE MIXES 
TEMPERATURE PROFILES. 



STABILIZES UNSTABLE 



SUBROUTINE MXXXT(TN.A) 
DIMENSION TN(20) f A (20) 
00 10 1=1 til 
IF(TN(I*t >.GE.TN(I 
IFUTN II 1-TN 1 1 *1 ») 
TAV = <TNCX*1)*TNI I) 



TN(I*1 )=T 
T\<I ) = TAV 
CONTINUE 
COMINJE 
TMAX CAMAX 



))G0 TO 1 
.LT. 0.0)00 
)/2. 



TO 1 



Ail 



(TN (1) ,TM2) , 
C,TN(6) ,TN(7 ) f TN<3) ,TN 
C(9),TN(10 ) ,TN(ll)t TN(12) ) 

IF(TN(12).LT.TMAX) GO TO lOu 

RETJRN 

END 



TN«3),TNim ,TN( 5) 



HO 



PLOTTER SYM CREATED ON 1 <! AJG 80 AT 12:56:46 

1 PARAMETER N-I 4 ,NN-12, NTlME -12, ND-110 

? 01 KENS I ON IBuFUOOGl 

} DIMENSION TEMP (501 ,DEE P( 5 J , TEMPS! ND I , DEEPSIND I ,QP (NM ,T2 INNI 

c S'iEffcfJSH L INI • AV r« , iCB(N»,2(NJ,XK2lNJ,7E0IND),THF(NDI t TSUINDl 

5 DIMENSION RGJ(N1 ,lMNl t DMCN».T2(NJ.ACNl,2EDCN0J 

6 OIHENSION Ell50),E<;(5Ql,E3(50»,E4l50J,E6l50!,E5J50J t 

7 CE7(50),E0(5O) ' 

8 CHARACTER*6 MONTHS IN) 

9 CHAR ACT ER*o IB CD 
JO H-l 

11 L:Q 

1? C 

13 C 

\1 C READ JRUN=1 IF YOU DESIRE PLOTS FOR MEASURED OAT* 

15 C READ JHUNr2 IF YOU DO NOT. 

1£ C NOIE : IF PLOTS FOR SEVEN STATIONS ANE N01 

17 C AVAILABLE, LINES 35 TO 16 MUST bE MODIFIED 

18 READ 100.JRUN.JYEAR 

19 100 FORMAT (I 

20 I COUNTED 

21 X2D=G. 

22 5 JO-0 

23 CALL PLOTS (IBUF , 1000, 1 1 1 

24 CALL PLOT(0. 0.7. 0,-31 

25 DO 1 Irl NTIME 

4* CALL REA0ER|T,AV.CtJ,2,A l XKZ,R0y t TN t OM i T2 t MONTHS t T2 f QP t 

H CCP.SIGMA.RS.RM.RS.Re.RV.Re-RQ.Rld.iPzlFREVLL.ROUCplDT! 

|8 «K20 t TE,N0AYS,f»a2 f Tl2 f Fi f F2,F3 f Fh i F41,F5 i F6.F7,FB,T6,TD2 f 

29 CNDAYS1 *TIM L, T 1ME2 , T IME3, I YL AR , MJ, XK . TDO,\;i ' 

30 I COUNT -J COUNT *1 

31 IFllC0UNT.iiT.96)60 TO 333 
3? IF (JRJN.EQ.2JG0 TO 200 

33 READI5 ,e > ( DEEP ( INK j, TF MP( I NK » , INK- 1, NS TOP ) 

3« DO 15 KL=1 ,50 

\l c READ L l5 I a^uEEP(KL) l El(KL) t L2(KLl,E3(KL),E4<KLl,ESlKLl t 

37 ^ADlbleiAEl.BEl.CEl.DEl.ELl.FEl.CEl.HLl.OEl 

Jo DEEP(KL1-A£1 

39 EllKLirBtl 

40 E2(KLI=CF1 

41 E3(KL)=DE1 
4? E4(KL)=EEI 

43 E5(KLJrFEl 

44 E6(KL1=6E1 

45 E7IKLJ=HE1 



16 EO(KL)=OEl 

17 IF(E3(KL1.CQ.Q.0IG0 TO 16 

18 IFIDEEPIKL J.tQ.I-1 .IIGO TO 16 

«*9 TEMP iKL)=1EliKL)«E2(KLI«£3(KL)«EllKL>*L5(KL>«Eb(KL)« 

50 CE7 (KL) I/.HHKL1 

51 TEMP IKL)=E3(KL) - 
5? 15 CONTINUE 

53 200 CONTINUE 

51 16 NST0P=KL-1 

55 IF URUN.EG.21G0 TO 201 

56 00 222 JIJ^l.50 

57 IFCDEEPIKL ).EG . ( -1 . ) )G0 TO 223 

58 READ(5,8)A£l,BEl ,CE1 ,t)El ,EE1, FE1 t GEI, HEl ,0E1 

59 IFIAE1.EQ, 1-iJ !G0 TO 223 

60 222 CONTINUE 

61 223 CONTINUE 

62 201 CONTINUE 

63 C0NS2=1./0.3Q18 

61 IF(JRUN.EQ.2)G0 TO 202 

65 00 9 INK = 1 ,NSTOP 

t6 DECP (INK )=L0NS2*DE£P( INK) 

67 9 DEEPIINK 1=150. -OEEPI INK) 

68 DEEP (NST0P*1)=Q.0 

69 DEEP (NST0P*2)=2(NN)/1.5 

70 TEMP (NST0P*11=0.0 

£ 71 TEMP (NST0P*2)=30. 0/1. S 

72 202 CONTINUE 
73 
71 
75 
76 
77 
78 

79 TEMPS (L)=T£MP (1) 

HO TCQIL)=TE 

£1 THF ILJ.= (Tl7l*T<eil /2. 

6 2 I BCD=MONTHSCJO) 

63 Z(NN*1)=0.0 

61 Z (NN*2 )=Z (NN1/1.5 

85 T(NN*1)=0.U 

66 T INN*2 )=30./1.5 

67 CALL AXIS 10.0 ,0.0 t 6HTEMP. ICIt-8,1. ^,0.0,1113) .Tlllll 

68 CALL AXIS(G.0,Q.0,9HDEPTH(FT),9, I. 5, 90. 0,7(13), 2(1111 
89 CALL FLINE(T,?,-NN, 1,0,0) • ' 

SO IF CICOUNT.GT. 96)60 TO 111 

91 IFfJftUN.E0.2IG0 TO 203 

9? CALL DASHL (TEMP ,0£ EP, NSTOP, 1) 

93 203 CONTINUE 

91 111 CALL SYNB0LI1. 0,0. 5,0. 11, 1BCD, 0.0,6) 

«5 CALL PLOT (2.25 ,0.0,-3) 



6 


FORMAT ( ) 


333 


JG=J0*1 




L=L*1 




TSJ(L)=T (12 




XZ0=XZD*30. 




ZEO(L)=XZO 



96 IFCJ0.EQ.4.0R. J0.Ew.81G0 TO 3 

97 GO TO 1 

98 3 CALL PLOT 1-9. 0,-2. 2S, -3 J 

99 1 CONTINUE 

100 CALL PLOT (-2.25,0.0,-3 1 

101 CALL SYHBOLI-6.75, b.75,.14,4lHTEHPERATLRE PROFILES FOR LAKE KEOUEE 
10? C ,0.0, till 

103 Pl-JYEAR 

104 MVrjYEAR 

105 CALL NUMBERI999. ,999. ,0.14 , PI , D.U.U1 

106 CALL SYMB0L(-6.75,b.S,0.1 t 54H(DEPTH IS MEASURED FROH THE DEEPEST P 
1C7 COINT OF THE LA KEI , 0.0 .54 1 

108 CALL PLOT IB. 0,-9.25,-31 

109 PRINT 2, MY 

110 2 FORMAT (IX, 'THE PLOTS F0R* t I5 t » ARE COMPLETE 1 ! 

111 IF CM. EO. 91 GO TO 6 

112 H=M»1 

113 JVEAR=JYEAR*1 

114 GO TO 5 

115 6 CALL PLOT (6. 0,0. 0,-31' 

116 00 13 1=1,96 

117 13 0EEPS(I1=2ED(I1 

118 DEEPS (97 1-0.0 

119 0EEPS(981=324Q. 0/9.0 

120 TSU(1091=Q.O 

121 TSU(H0)=35./5. 

122 TEQ(l09irn.0 

123 TEG(1101=35./5. 
121 THF(109)ra.G 

125 THF (1101=35. /5. 

126 TEMPS (97 1=0.0 

127 TEMPS (98 1=35. /5. 

128 ZEb(109l=0.G 

129 ZED(1101=3240./9. 
120 CALL PLOTIG. 0,2. 0,-31 

131 CALL AXIS (0.0,0.0, 12HTIME I N DAYS, -12, 9.0,0.0 , ZED (109 1 ,ZE0 (1 101 1 

132 CALL AXIS (0.0,0.0, 16HTEHPERATURES (CI , 16, 5.0, 9u. , TSU (109 1 ,TSU 

133 CU1G11 - 

134 CALL FLINE|7E0,TSJ,-108,1.2,111 

135 CALL FLINEIZED,TE0,l08,l t 2,2) 

136 CA!.L FLINE(ZED, THF, -106, 1,2,01 

137 IF4JRUN.EQ.21G0 FO 204 

138 CALL UASHL (DEEPS ,TLMPS ,96 , 11 

139 . 204 CONTINUE 



140 c 

111 c 

142 C 

1 43 C CHANGE TITLES TO SUIT NEEDS (4 LINES* 

144 C 

145 C 

146 CALL SYMBOLIC 0,6.0,0. 14. 

147 C46HSTRATIF1CA1I0N CYCLE FOR LAKE KEOkEE 1971-1979,0.0,461 

148 CALL SYMBOL (0.0,0. 10, 0.1O,67H 1971 1972 1973 1974 

149 C 1975 1976 1977 1978 1979,0.0,87) 

150 URITEI6.71 

151 7 FORMAT UK . 'ALL PLOTS ARE NOW COMPLETE •,//, * NORMAL JOB EXIT'l 

152 CALL PLOT 115.0,0.0,-31 

153 STOP 

154 EriD 






Ul 



READER SYM CREATED GN 12 AUG 80 AT 13:21:15 

1 C 

2 C 

3 C THIS SUBROUTINE READS THE MAGNETIC TAPE 

4 C CONTAINING THE COMPUTED RESULTS. 

5 C 

6 C 

7 SUBROUTINE READER I T.aV.CB t 2,A, XK2, HOW, Tn, DM, TZ. MONTHS ,T2,QP, 

8 CCP ( S1GMA,R3,R<< ,R5 t R6,R7,R8.R9.RlD L GP2.FREVEL.R0VCP.DT. 

9 CXK20 I TE,NDAYS 1 TN12,T12,Fl 3 F2,f3,F3l,F<»l t FS.F6 l F7,F8 1 T0.TD2, 

10 CNDAYSl ,TlME,TlME2, TIME 3,1 YEAR, MJ, XK, TOO, J,NCASE .SF.EDEPT ,VOL I 

11 DIMENSION TI201 ,AV»20I .CBJ20) ,Zf 20) ,A ( 20 I , XK2I 20 >. ' 
1? CROb<?0 ),TNf20> ,DM< 201, 12*2 01, I 2f 2Q) ,QP 1121 

13 CHAR ACT ER*6 MONTHS 1 121 

14 1 CONTINUE 

15 READ (6,END=11 C T< Ul , IJ=1 , 121 , I A V<U> ,U=1 ,12 I. 

16 CtCblU),lJ = l l 12) .(2CIJ1.IJ-1, 121. <A(IJ),IJri.l2), 

17 CIXK2 (IJ) ,U-I ,12), (R>0UllJl,lJ=1.12j,lTNffIJI,lj=l v |2lj , 

18 C(DH(lj),ljri,i2) ,C T2C I Jl , U-l, 121, CMGNTHSf Ul , U=l,12 I , I 

19 CIT2IUI,IJ-1,121 , ' 

20 CTQPfU l,U=1.121 , 

21 CCP,SIGHA,R3.R<i ,R5,R6 t R7.R8,R9«RlU,CP2,FREVEL,RCUCP v DT 9 

22 CXK20 t TE,NDAVS 1 TM2 1 Ti2.Fl,F2,F3,F3l,Fll t F5,F6 1 F7 1 Fa,T0,T02, 

23 CNDAYSl ,TlMt,TIME2, TIME3, I YEAR ,MJ, XK, T00,J,NCASE ,SF .EOEPT ,¥0L 
2«* RETURN 

25 ENO 



O) 



SHOOTH SYH CREATED ON 12 AUG 80 AT 1H:31:30 

1 C 

2 C '.HIS SUBROUTINE CORRECTS THE EDDY UlFFbSlVITY 

3 C IF VARIABLE TIME STEP IS REQUIRED, •DTl« SHOULD 
1 C BE CHANCED TO *DT * IN THE CALLING PROGRAM. 

5 C 

6 SUBROUTINE SHOOTHI XKZ,XK2U ,XK2L,NDAYS 1 , TN12, Tl2tT, Ofl ,D2l 

7 DIMENSION XK2 (20 1 , Tl 20 I 

8 DO 93 1=1,12 

9 IFIXKZUl.CT.XKZUl XKZdjrxKZU 

10 IF (XKZ (Il.Ll.XK2D XKZII1=XKZL 

11 93 CONTINUE 

12 Nti.ro 

13 00 96 1=2,12 

11 IF (XKZ(I).EG.XKZLl NEH=I 

15 96 CONTINUE 

16 IFINEW.E0.01G0 T0?7 

17 00 55 1=1 .NEW 

18 XK2lI)=XKZL 

19 55 CONTINUE 

20 77 CONTINUE 

21 IF(N0AYS1.LE.60.0R.NDAYS1.GT.3G0) GO TO 29 

22 IFUN12.GE.TI21 GO TO 19 

23 IF UN12.LT. 1121 GO TO 39 

21 19 XMIN^AMINI <XKZ (1) , XKZ 1 2 1 , XK2< 3 1 ,XKZ( 1 1 ,XK2 (5l , XKZ tfc I ,XKZ (7 1 ,XKZ 

25 1(8).XKZ(9>,XKZ(10I,XKZ(111,XKZ(12I> 

26 DO 82 1=1 ,12 

27 IF(XKZlIl.EQ.XMlN) GO TO SI 

28 82 CONTINUE 

29 GO TO 29 
JD 81 IMIN=I 

31 00 70 1=1 ,1HIN 

32 XKZ(I|=XK2(IMIN) 

33 70 CONTINJE 
31 GO TO 29 

35 39 XMAX-AMAX1 |XK2 1 1 ) , XKZ I 2 1 ,XKZt 3 I ,XKZl 1 1 ,XKZ» 51 f XKZ 161 ,XK2 C7 1 ,XK2 

36 l(8),XKZ(9>,XKZ(10l ,XKZ ( 1 1 1 ,XK2 ( 1 21 1 

37 00 62 1=1 .12 

38 IF tXKZtll.EG.XHAXl ISO TO 61 

39 62 CONTINUF 
«*0 GO TO 29 
11 61 IHAX=I 

**2 DO 50 J=IHAX,12 

13 XKZ(II=XKZ(lMAX) 

»*1 50 COfiTlNjE 

15 29 CONTINUE 

16 200 XMAX-AMAX1 (XKZ 1 1 1 , XKZ I 2 1 , XiiZI 3),XK2(<<I , XKZ < 51 , XKZ (6 ) ,XKZ 17 ) ,XKZ 

17 lffcl,XKZ(9l,XKZI101 ,XKZ(111 ,XKZ(12)) 

18 DT1=(U.1*DZ**21/XHAX 

19 RE URN 
50 END 



wp*r 



STORE SYM CREATEL ON 12 AJG 60 AT 13;19:<i7 

1 C 

2 C 

3 Z THIS SUBROUTINE STORES THE COMPUTED RESULTS ON 
1 C MAGNETIC TAPE. 

5 C 

6 C 

I SUBROUTINE STOPECT.AV, CB, 2, A, XK2,R0k,TN f OH, T2 f K0NTHS,T2,QP, 

* LX RZO.lt, NO AYS , IN 1 Z ,T1Z t F 1 , F 2. F 3, F 31, F «♦ 1 , FS.Ffc, F7 , F£ , TO ,TD2 . 

12 

13 COP U 12 I 

1*» CHARACTER*^ MONTHS (12) 

}5 U'RITE C61 ITIIJl ,Iu:l,12S,|AVIlJl,IJ=l,121, 

J 6 , Sf^il-J^flW^illflt^jllJl.ljri.lijjAllJl.lJri,! 

18 
19 

2n ccQpiiji;u=i:i2i , 

\\ ^S^^S I ^V^^i? l, *?.?i , ^!?l 7 * Ra A5 9 i5 1 yi.1: p 2 l FREVLL l RpWCP f 0I 

23 




^ 2H END FILE 8 

25 RETURN 



CJCBIIJI.IJ=1.12» tf 2UJI.IJ= 1,1 21.1 AflJ).lj£l t 121. 
CIXK2lIjJ t IJrI,i2l f (R0WllJl,ljri,X2i,cilJlJJ t Ij = i;i2j, 
C«DKIIJ),IJ=1,12) ,( 121 IJ1, IJ=1, 121, IMONTHSdJI.ljr 1,121, 
C( 1 2 1 1 J 1 , 1 »J— 1 ,121, 
CCQP(IJ1,IJ=1,12] - 

^S^»,S I ^* , ^ :j »? l, * R ?i ,<6 ! n7 * Ra A R9 i ftltJ i < - P2 F Fft E: v tL l R0WCP f 0T 1 
CXK20 ,TE, N0AYS r fN12,Tl2,Fl,F2,F3,F31,F&I,F 5 t Ffc,f7,FS ,10,702, 
CNDAYSl,TlME,TlME2,lIHE3,IYEAR,HJ,XK,TOL,J,NCASt,SF,£OEPT,VOI 



VOL 

ULt H 
_ - JRN 
26 END 



I 

2 
3 

4 
5 
6 

7 

a 

9 

13 
11 



YEARS SYM CREATED ON 12 AUG 80 AT 1J:1U:C3 

C 

C THIS SUBROUTINE PRINTS THE YEAR TITLE. 

C 

SUBROUTINE YEARS ( S ELTEIM, OQPP, I YEA H ) 

PRINT 99.IYEAR 
99 FORMAT 159X ,17 <•*• I ,/, 5 9X,'**, 1SX, ***,/, 59X, 

C*** ,2X ,*YEAR - NI^tEX, •*•,/, 59X, •*», 15X f « »• 
C./.59X »17( '*') > 

RETURN 

ENO 



48 



